Conductance of p-n-p graphene structures with 'air-bridge' top gates 
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We have fabricated graphene devices with a top gate separated from the graphene layer by an air 
gap — a design which does not decrease the mobility of charge carriers under the gate. This gate is 
used to realise p-n-p structures where the conducting properties of chiral carriers are studied. The 
band profile of the structures is calculated taking into account the specifics of the graphene density 
of states and is used to find the resistance of the p-n junctions expected for chiral carriers. We show 
that ballistic p-n junctions have larger resistance than diffusive ones. This is caused by suppressed 
transmission of chiral carriers at angles away from the normal to the junction. 



Graphene, a monolayer of carbon atoms, is a new two- 
dimensional system [1 with unusual properties: the gap- 
less energy spectrum of electrons and holes is linear, and 
these carriers are chiral. It is the chirality that suppresses 
the backscattering of carriers ^ and allows them to pen- 
etrate through potential barriers without reflection (the 
Klein paradox [5]). The theory of ballistic graphene p-n 
junctions jl] predicts that the carriers propagate through 
the barrier without scattering if their angle of incidence is 
0° (with respect to the normal) and are partially reflected 
at other angles. This angular selectivity of the carriers 
determines the resistance of a ballistic p-n junction. 

To realise ballistic transport through a graphene p-n 
junction, the mean free path of carriers / has to be larger 
than the characteristic length of the junction 2t. A num- 
ber of intriguing phenomena, such as small-field positive 
magnetoresistance j4j and oscillating transmission prob- 
ability [3], can be observed in a ballistic p-n-p structure 
where the total length (two junctions plus the n- region) is 
smaller than /. Ballistic p-n and p-n-p structures can also 
be important for a number of potential applications: for 
example, graphene lenses [5] and filter circuits [3] where 
the electrical current of chiral carriers is focused on or 
directed to a desired contact. Therefore, exploring the 
ways of producing high-mobility p-n-p graphene struc- 
tures with large / and understanding the mechanisms of 
carrier propagation through them is an important task. 

Graphene flakes obtained by mechanical exfoliation 
[1] are conventionally deposited on a Si02/Si substrate 
where the conducting n-Si layer forms a (back) gate act- 
ing on the whole flake. Inversion of the type of carrier 
in a part of graphene flake has been recently achieved by 
using an additional (top) gate which is positioned above 
the graphene layer [6j [7l [8l |9] . In the p-n junctions fab- 
ricated so far, the top gate rests on an insulating layer 
that can decrease the mobility of graphene carriers un- 
der it [B] and hinder the realisation of ballistic p-n and 
p-n-p structures. We avoid this problem by fabricating 
p-n-p graphene structures using suspended 'air-bridge' 
top gates. (Similar designs were used earlier in semicon- 
ductor nanostructures [10l|TT].) The mobility of carriers 
under the top gate in such structures is the same as in 



the rest of the graphene layer. 

The theory of ballistic graphene p-n junctions [4j con- 
siders a 'smooth' junction, 2k pt ^ 1, where kp is the 
Fermi wave vector of the particles and 2t is the tunneling 
distance of the carriers. The transmission probability of 
such a junction as a function of the angle of incidence 9 
(with respect to the normal to the junction) is given by 
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Here, w^=10^ms^^ is the Fermi velocity of carriers and 
F/e is the electric field in the barrier which is assumed 
to be constant over the tunneling distance so that the 
potential ip{x) — Fx/e. The term ky is the component 
of the wavevector parallel to the junction, ky = nn/W, 
n — 1,2, . . ., where W is the width of the barrier. There- 
fore, the transmission of chiral carriers is restricted to 
the angles 9 < 9c — {F/Trhvpkp)^^''^ or, equivalently, the 
transverse momentum values ky < {F/ttIIvf)^''^. The 
conductance of a ballistic junction is then [4] 
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For the calculation of the resistance of a graphene p- 
n junction it is important to know accurately the elec- 
tric field at the p-n boundary, which can be much larger 
than that in the rest of the barrier [12] . This is caused 
by poor screening near the point of electro-neutrality 
(the Dirac point) where the density of carrier states de- 
creases to zero. Therefore, to find the expected resis- 
tance of ballistic p-n junctions in our samples, we have 
computed the potential profile of the experimental p-n- 
p structures. This is done for different combinations of 
the back and top gate voltages, taking into account the 
density of states in graphene which changes linearly with 
energy. 

To fabricate ballistic p-n junctions, one positions the 
top gate close to graphene, which increases the field 
F and decreases the tunneling distance t{ky,F) — 
hvpky/F. In addition, the Fermi wavelength increases 



when the electron diffusively approaches the junction, 
up to the distance ~ 1/2 from its centre when it passes 
through it ballistically. As a result, the length of the 
junction 2i can become comparable to the wavelength, 
and the assumption of |1] of a smooth barrier not ap- 
plicable. To examine the validity of this approximation 
for our structures, we have also found the transmission 
probability w{9) by solving the Dirac equation numeri- 
cally [13], using the calculated potential profile <^{x). 

This calculated potential ^p{x) allows us to determine 
the resistance of the p-n-p structure expected for diffu- 
sive propagation of carriers (without taking into account 
their chirality) and to compare the result with exper- 
iment. The required resistivity p{^p) of the sample at 
different Fermi energies is found from the resistance as 
a function of back-gate voltage Vbg (at top-gate voltage 
Vta = 0). 

Here, we present the results for three samples with dif- 
ferent mobilities. The sample with the smallest mean 
free path has a resistance in the p-n-p regime which is 
in agreement with the diffusive model, while the higher- 
mobility samples (with ballistic transport through the 
p-n junctions) show an enhanced value of the resistance, 
which agrees with the result we expect for the chiral car- 
riers. 
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FIG. 1: (a) Three stages of the air-bridge fabrication: elec- 
tron beam lithography with two exposure doses; development, 
and deposition of the metal film, (b) A false-colour SEM im- 
age of a graphene flake with a metal air-bridge gate, tilted at 
45°. (c) Grey-scale of the positions of the maxima in R as 
a function of carrier density and magnetic field, with dotted 
lines corresponding to the shifts of the lowest Landau levels 
expected for single-layer graphene. (d) Positions of the resis- 
tance maxima at different Vbg, Vtg, with different regimes of 
the device operation indicated, sample SI. 



To fabricate the p-n-p structures, we have chosen 
graphene flakes of rectangular geometry on Si02/Si sub- 
strates with a 300 nm oxide layer. The samples have the 
following dimensions, in fim: L=5, VF=0.24 (sample SI); 
i=4.3, W=0.6 (sample S2) and L=1.45, W=0.15 (sam- 
ple S3). The mobility of these samples outside the region 
of electro-neutrality (at a carrier density of 3 x 10"'^^ cm~^) 
is 13, 11 and 6 x lO'^ cm^V~^s~^, respectively. The proce- 
dure of the top gate fabrication is illustrated in Fig. Ilja) . 
Two layers of PMMA with different molecular weights are 
spun on the flake: a soft resist (495K) on top of a hard 
resist (950K). They were then patterned using 10 kV e- 
beam lithography (to achieve larger undercut in the top 
PMMA layer). Two different exposure doses were used 
in the areas of the span and pillars of the bridge, while 
the area outside the bridge was not exposed. The dose 
in the span is just enough to expose the soft resist, but 
too small to affect the underlying hard layer. Both lay- 
ers are exposed at a larger dose in the areas of the pillars 
(and contacts). The structures are then developed and 
covered with 5/250 nm of Cr/Au. The 'lift off' removes 
PMMA leaving the bridge with a span up to 2 /im sup- 
ported by two pillars. Figure W[h) shows an SEM image 
of sample S2 with a bridge top gate and two Ohmic con- 
tacts. The mean free path in our samples / « 45— 100 nm 
and the distance between the top gate and the flakes is 
130-210 nm. In the attempt to produce a p-n-p struc- 
ture with ballistic properties, the top gate is made short 
in the direction of the current flow: 100-170 nm. 

Two-terminal measurements of i?(Vbg) in quantising 
magnetic fields have confirmed that we are dealing with 
single-layer graphene. The grey-scale plot in Fig. IIJc) 
shows the shift with magnetic field of two resistance 
peaks corresponding to the first electron and hole Lan- 
dau levels. In this experiment, the carrier concentra- 
tion was varied by the back gate voltage: n/Vbg — 
7.2 X 10^°cm~^V~^ (the relation from the known capac- 
itance of the structure with a 300 nm Si02 layer). A 
positive (negative) Vbg induces electrons (holes) in the 
graphene layer, and the 0th Landau level (the bright 
vertical line) corresponds to the Dirac point. The dot- 
ted lines show the position of the 1st Landau level of 
electrons and holes in accordance with the filling factor 
1/ = 4eBn//i = ±1 [Ml HI]. 

Figure [nd) shows that, by changing the combination 
of the voltages on the back gate and top gate, one can 
get different regimes of the operation of the device, with 
both p-n-p and n-p-n modes available. Different regions 
are separated by the resistance peaks corresponding to 
the Dirac points under the top gate (steep line) and back 
gate (almost horizontal line). The slope AVbg/dVtg of 
the steep line (from 0.24 to 0.4 in our samples) gives the 
efficiency of top-gate control with respect to that of the 
back gate. 

Figure [2F a) shows R{Vbg) at Vtg = 0, where the peak 
corresponds to the Dirac point. Figure l2{b-d) shows how 
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FIG. 2: (a) Resistivity of the three samples as a function 
of the back-gate voltage, at Vtg = 0. Points indicate the 
fixed values of Vtg where the top-gate voltage was swept to 
produce p-n-p junctions, (b) The resistance of sample SI 
as a function of top-gate voltage at different Vbg. (c,d) The 
resistance as a function of top-gate voltage at different Vbg of 
samples S2 and S3, respectively. Points show the results of 
the calculations of the expected resistance assuming diffusive 
transport of carriers. (Dashed lines in b,c are guides to the 
eye.) 



the resistance of the samples changes (for fixed Vtg values 
shown by points in Fig.l2ja)) when a top-gate voltage is 
applied. When the main part of the sample outside the 
top gate is p-type, applying a negative Vtg decreases the 
resistance due to the increase of the hole density under 
the top gate. Applying a positive Vtg increases the re- 
sistance of the samples, first because of the depletion 
of electrons under the top gate and then because of in- 
version of the sign of carriers and formation of a p-n-p 
structure. 

At the onset of the formation of the p-n junction, the 
resistance shows reproducible oscillations as a function of 
Vtg. They survive at high temperatures, which strongly 
suggests that they can be due to the oscillations of the 
transmission coefficient caused by interference of chiral 
carriers [3] within a ballistic p-n-p structure. However, 
this effect has to be separated from mesoscopic fluctua- 
tions of resistance |16j that can be enhanced in the Dirac 
points of the p-n junctions. (We will discuss the sep- 
aration of these effects elsewhere [l2]-) Here, we con- 
sider the average values of the resistance and present the 
results obtained in a higher temperature range (r=50- 
77 K) where the effect of the mesoscopic fluctuations is 
less important. 
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FIG. 3; (a) Geometry of the top-gated structure used in the 
calculations of the potential profile along the flake, sample 
SI. (b,c,d) Potential proflle of the top-gated samples SI, S2 
and S3 along the barrier at different Vtg for a fixed Vbg. The 
curves correspond to the position of the Dirac point and zero 
is the Fermi level. The bars show the mean free path /. 



The band-structure profile along the p-n-p structure 
at different Vtg is calculated numerically for the geome- 
try shown in Fig. p[a), by solving the two-dimensional 
Laplace equation with the potentials of the two gates 
as boundary conditions. The presence of the flake is 
included as an additional boundary condition on the 
jump of the normal component of the displacement 
field at the flake: AZ?„ = en(x), where n(x) is car- 
rier concentration along the sample. Taking the Fermi 
level as zero and using the linear energy dependence of 
the density of states in graphene, i^(-E) = 2E / ttIt? v^p , 
one gets the relation between the carrier concentration 
along the junction n{x) and the electrostatic poten- 
tial: n{x) = sgTi{ip)e'^ip'^{x)/-K{hvpY , where sgn((p) re- 
flects the fact that carriers can be both electrons and 
holes depending on the position of the Dirac point with 
respect to the Fermi level. Examples of such calculations 
for the three samples are shown in Fig. p[b-d) , where in- 
deed one can see a rapid increase of the electric field at 
the boundaries between p- and n-regions considered in 
[12]: i^ = (0.8 - 2.4) X 10^ eV/m. 

We use the 'cahbration' curves R{Vbg) in Fig. Wia) to 
find the resistivity of the flake at different Fermi ener- 
gies (different electrostatic potentials if the Fermi level is 
taken as zero). We used the relation between Vbg and the 
electrostatic potential ip known from the capacitance be- 
tween the flake and the back gate: (/^(mV) = Sl-^/VbgTV)- 
With the calculated distribution of the potential along 
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FIG. 4: (a) Momentum of the electron approaching the junc- 
tion at an angle 9. (b) The gap in the spectrum E{kx) (high- 
lighted) at ^ 7^ 0. (c) Band-structure profile along the length 
of the p-n-p structure. The value of the gap determines the 
tunneling length 2t(6). 



the structure ^p{x), the integration of the resistivity p{Lp) 
gives the resistance expected for diffusive propagation of 
carriers: R = {^/W) J p{x) dx. 

In the range of Vtg corresponding to accumulation 
(negative Vtg) and depletion (small positive Vtg) under 
the top gate, the resistance is well-described by the dif- 
fusive model, Fig. [2Fb-d). One adjustable parameter, the 
distance h between the top gate and the graphene flake, 
was used in plotting the calculated values: h =140, 210 
and 130 nm for samples SI, S2 and S3, respectively. The 
obtained values are close to those expected from the fab- 
rication process and agree with observed efficiency of the 
top gate. Fig. |2jd). With larger positive Vtg and for- 
mation of the p-n-p structures, samples SI and S2 show 
significantly larger values of the resistance than expected 
from the diffusive model: Ai? ~4 and 2kOhm, respec- 
tively. However, the narrowest sample S3 with the lowest 
mobility shows agreement with the diffusive model in the 
whole range of Vtg , Fig. WiA) ■ 

To explain these observations, we find the character- 
istic thickness of the p-n junctions in the three sam- 
ples and compare it with the mean free path I. Ac- 
cording to [4], the reason for the enhanced resistance 
of a junction is the decrease of the transmission when 
the electron approaches the junction at an angle 9 ^ Q, 
Fig. l4{a). Conservation of the parallel component of the 
momentum ky produces a gap in the energy spectrum 
E{kx) for the motion across the junction, Fig.EFb). The 
distance 2t is then defined as the classically inaccessi- 
ble region which requires electrons to tunnel along it. 
Fig. Wc): t — hvphpsinO/F. The critical angle for car- 



rier transmission in the three samples varies in the range 
6'c = 20 — 30°, assuming the length of the ballistic p- 
n junction to be I and taking the fci?-value at a point 
X — —1/2 from the barrier. Fig. l4Fa). As the tunneling 
distance 2i depends on the angle of incidence, we take for 
a typical value of the barrier thickness 2t{29c) —40 nm in 
our samples. 

The mean free path / has been found using i?(T4g) of 
a uniform sample at Vtg = 0, Fig. l2ja), and the rela- 
tion (J — 2e^{kpl)/h. The value of I weakly depends on 
Vbg, and when extrapolated to the Dirac point (Vbg = 
for an undoped sample) gives / ~ 100, 75 and 45 nm, 
respectively, for samples SI, S2 and S3. Comparing the 
tunnelling length with the mean free path shows that the 
p-n junctions in SI and S2 are ballistic {I ^ 2t), while in 
S3 they are less ballistic {I ~ 2t). This can explain the 
agreement of the resistance of S3 with the result of the 
diffusive model in Fig. ^d) . 

To find the expected resistance Rpn of ballistic p- 
n junction in samples SI and S2 and compare it with 
the observed difference Ai? in Fig. |2] we first assume a 
smooth potential barrier, 2fci?i ^ 1, and by using the 
calculated value of electric field F we get the tunneling 
probability w{9) from Eq. [I] Equation |2] is then used to 
obtain the resistance of the ballistic p-n junction. We 
have found that using summation rather than integra- 
tion is more appropriate in our case, as samples SI and 
S2 have less than 12 modes (the narrowest sample S3 has 
only three modes). The value of the Fermi momentum 
kp in these calculations is taken at a distance 1/2 from 
the barrier using the values of the mean free path found 
above; however, the result for Rp^ hardly changes if the 
value of I is varied by two times either way. This is clear 
as the tunneling probability w{9) in Eq. n] depends only 
on ky which takes specific, quantised values ky = nn/W. 
The obtained values are Rpn =5 and 2 kOhm for samples 
SI (at Vbg = -9V, Vtg = 40 V) and S2 (at Vbg = -4V, 
Vtg ^ SOY). 

Taking into account the Fermi wavelength at the dis- 
tance 1/2 from the barrier, we see that 2fci?i ~2 for the 
three samples. To examine the applicability of a smooth- 
barrier approximation for this (not too large) value of 
2kpt, we have calculated w{9) directly using numerical 
methods fl3J and compared the result with that obtained 
from Eq. [l] It shows less than 5% difference from the 
value of Rpn calculated above and a significantly larger 
resistance than the one expected for a sharp, rectangular 
barrier where w{9) — cos^ 9 |1]. 

In experiment, it is not the resistance of an individ- 
ual ballistic p-n junction which is measured but the re- 
sistance of the whole p-n-p structure. It can be differ- 
ent depending on whether its middle, n-region is long or 
short compared with I (i.e., diffusive or ballistic). For a 
diffusive n-region with three independent contributions 
(two junctions and middle region) Rpnp > 2i?p„, while 



for a ballistic n-region, Rp 



Rpn HI- The resistance 



of a ballistic p-n-p structure should not increase with 
addition of another junction as the electrons approach- 
ing the second junction have already been selected by 
the first junction within the critical angle 9c- Therefore, 
they all will have high transmission probability 'w{9) go- 
ing through the second junction. 

Figure |2|b,c) shows clearly that the resistance of SI 
and S2 is larger than that expected in the diffusive model 
by Ai?, because of the ballistic transport of chiral carriers 
through two p-n junctions. To find their resistance, we 
assume that they are independent; that is, the n-region 
is diffusive. Then the observed difference Ai? — 2(i?p„ — 
i?^), where i?^^ is the resistance of the diffusive p-n 
junction on the length I which was taken into account 
in the diffusive-model calculation shown in Fig. [2j With 
the values I =100 and 75 nm, one finds that _R£j =2 and 
0.6kOhm for samples SI and S2, respectively. This gives 
the corresponding resistance of the ballistic p-n junction 
Rpn =4 and 1.6kOhm, which is close to the expected 
values of 5 and 2kOhm. (Even better agreement, within 
10%, is achieved if another quantisation rule for graphenc 
is used [18]: ky = 7r(n + l/2)/W, n = 0, 1,2, . . .) The 
assumption of the diffusive nature of the n-region at large 
Vtg is confirmed by Fig. l3| where the whole p-n-p region 
is seen to be larger than the mean free path. However, 
near the onset of the p-n junctions, at small Vtg, the p-n-p 
region is much shorter and can be fully ballistic. 

In conclusion, we have fabricated p-n-p and n-p-n 
graphene structures using non-perturbative 'air-bridge' 
top gates. The chiral nature of charge carriers in 
graphene has been directly demonstrated by detecting 
an increase of the resistance of p-n junctions caused by 
their selective effect on the propagation of chiral par- 
ticles. Our detailed analysis shows that individual p-n 
junctions are ballistic and that a ballistic p-n-p structure 
can be realised using this fabrication method. 
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